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AN  HF  COMMUNICATIONS  FREQUENCY-MANAGEMENT  PROCEDURE 
FOR  FORECASTING  THE  FREQUENCY  OF  OPTIMUM  TRANSMISSION 

I.  INTRODUCTION 

The  needs  Co  utilize  Che  crowded  tur  cnannels  efficiently  ana  Co  maintain  a 
hign  degree  of  communicaCions  channel  conneccivity  justify  a  substantial 
effort  co  accurately  predict  cne  values  of  those  ionospheric  parameters  from 
wnicn  optimum  frequencies  may  be  selected.  Oblique  incidence  sounders 
operating  over  all  necessary  communications  pachs  would  give  a  near-ideai 
amount  of  the  information  required  for  optimum  frequency  management.  In 
particular,  they  give  tne  Maximum  Observable  frequency  uHO f)  with  an  rms  error 
on  the  order  of  1  Mhz,  they  give  relative  amplitude  as  a  function  of 
frequency,  and  they  give  information  on  muitimouing.  But  tne  procurement  and 
operation  of  that  many  sounders  would  be  unrealistically  expensive;  sounders 
cannot  forecast  cuCure  values  of  Che  noF,  and  there  are  problems  associated 
with  their  use  in  denied  areas  and  in  covert  situations.  Vertical  incidence 
sounders  have  similar  advantages  and  disadvancages .  as  a  consequence, 
ionospheric  models  based  on  extensive  periods  of  ionospheric  measurements  are 
routinely  used,  ana  will  probably  continue  to  be  used  lor  a  long  time  to  come, 
to  predict  the  parameters  required  for  HF  frequency  selection.  One  of  the 
most  successful  of  these  is  the  MINIMoF  model  developed  by  the  Naval  Ocean 
systems  Center  (NObC).  This  prediction  algorithm  requires  only  very  modest 
computing  facilities  ana  is  thus  suitable  for  field  use.  Using  only  a  current 
(within  1*5  days;  value  of  the  solar  10.7  cm  ionizing  flux,  MINLMUF  is  capable 
of  predicting  the  MOF  over  any  desired  paths  with  an  rms  accuracy  of  3-5  Mhz. 
Other  models,  such  as  10NCAP,  are  slightly  more  accurate  but  require  much  more 
extensive  computing  facilities  and  are  thus  not  easily  field-usable.  None  of 
the  models,  however,  gives  forecasts  with  an  accuracy  approaching  that  of  the 
sounder  measurements.  ihis  leads  u8  to  attempt  an  improvement  in  forecasting 
accuracy  by  combining  sounder  measurements  with  predictive  modelling.  With 
sounder  data  from  a  small  number  of  sounders  used  as  input,  the  models  can 
Chen  produce  predictions  for  a  much  larger  number  of  communications  paths. 

The  goal  is  co  produce  predictions  with  an  accuracy  approaching  that  of 
sounder  measurements  and  with  convenience  approaching  that  of  model 
calculations.  This  is  the  basis  of  the  NRL  "updating"  technique.  It  is  the 
purpose  of  this  technical  note  to  investigate  the  foundations  of  the  Updating 
tecnnique  and  to  derive  one  of  its  numerous  possible  applications:  a 
procedure  for  using  the  MOF  measured  on  a  sounded  path  to  forecast  the  highest 
frequency  on  an  unsounded  path  wnicn,  with  a  specified  degree  of  probability, 
will  not  exceed  the  MOF  on  that  path.  This  frequency  is  called  the  frequency 
of  optimum  transmission,  or  FuT. 

To  illustrate  tnese  techniques,  data  taken  from  tne  bolid  Shield 
experiment,  a  test  of  HF  communication  procedures  in  the  Eastern  United 
States,  were  used.  These  data  were  selected  because  they  covered  a  relatively 
long  time  interval  (17  days)  and  were  available  in  a  readily  usable  form. 

Work  on  Updating  procedures  of  NKL  has  used  data  from  oblique-incidence 
sounders.  Only  1-hop  propagation  modes  have  been  considered,  so  that  tne  MUF 
values  obtained  are  characteristic  ot  the  ionosphere  at  the  midpoints  of  the 
paths,  bimilarly,  the  correlation  of  properties  on  the  two  sounded  paths 
actually  measures  the  correlation  or  ionospheric  properties  ac  the  midpoints 
of  the  two  paths.  These  conditions  differ  from  the  conditions  present  in 
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ocher  investigations  referenced  in  this  reporc,  in  which  only 
vertical-incidence  sounder  data  were  used. 

The  correlation  of  ionospheric  properties  on  two  patns  separated  in  time 
and/or  space  must  be  made  through  a  model  which  incorporates  the  temporal  and 
spatial  MUr  uependence  ana  wnicn  can  oe  adjusted  by  a  time-anu-space  invariant 
parameter.  In  the  present  worn  tne  MliilMUF  model  is  used  exclusively  because 
of  its  simplicity  anu  relative  accuracy.  The  moael  is  aajusted  by  a  single 
parameter,  the  10.7  cm  solar  tlux. 

Section  II  of  this  report  contains  a  definition  and  examples  of  the 
correlation  coefficient  connecting  the  modelling  parameters  measured  on  a  pair 
or  paths,  and  discusses  now  tne  correlation  coefricient  depends  on  the  time  of 
day. 

In  Section  III  three  examples  of  torecasting  are  discussed.  In  the  first 
case,  MINIMUF  is  used  alone  tc  forecast  the  MllF.  In  the  other  cases  the 
M1MIMUF  prediction  of  MoF  on  a  specified  path  is  modilied  by  using  .nOF  values 
measured,  respectively,  on  the  same  or  a  different  path.  The  quality  of  the 
forecasts  is  measured  by  comparing  the  rms  differences  of  the  rorecast  and 
measured  values. 

In  Section  IV  it  is  snovn  how  to  forecast,  for  an  unsounded  path,  the 
Frequency  of  Optimum  Transmission  -  that  is,  the  highest  trequency  which  is 
assured,  witn  a  specified  probability,  to  De  below  tne  MOF  on  that  path,  when 
the  MOF  on  a  sounded  path  and  the  correlation  of  the  modelling  parameter  for 
the  two  paths  are  known. 

II.  MOF  CORKtLATlOiN 

Correlation  Parameter 

While  the  Mu*  is  consiaered  the  single  most  userul  ionospheric  property 
for  frequency  management  purposes,  one  cannot  directly  correlate  MOFs  at - 
different  positions  on  the  eartn  because  of  the  MOF's  well-known  dependence  on 
latitude  ana  time  of  day  (Rush  1976).  Instead,  we  adopt  a  MOF-prediction 
model  which  incorporates  the  effects  of  geographical  positon  and  time  of  day, 
and  look  for  some  parameter  which  can  be  adjusted  to  make  tne  predicted  MOF 
equal  the  real  li.e.,  measured)  MuF.  Tne  MINIMUF  model  (nose  et  al  i978a,b) 
was  chosen  because  it  is  easy  to  use,  reasonably  accurate,  and  has  enjoyed 
considerable  success  in  frequency-management  operations.  MINIMUF  computes  the 
MuFs  by  the  equation 

MOF  -  (1  ♦  K/250)  M  ♦  A^  Icos  Xeff7 

where  k  ■  sunspot  number,  Aq  and  are  constants,  Xefi  i®  c^e  effective 
zenith  angle,  and  M  is  a  function  of  pathlength  ana  ionospheric  layer 
structure.  The  MINIMUF  sunspot  numoer  R  is  obtained  from  the  values  of  solar 
10.7  cm  flux  F0  broadcast  by  WWV  according  to  the  formula 

R  -  56i.d  j[10. 52996  -  0.U00J56  (6J.75  -  F0)J  -  0.72b| 

The  sunspot  number  so  calculated  is  nearly  directly  proprotional  to  the 
lu.7  cm  flux  F0.  following  recent  custom,  we  use  tne  10.7  cm  flux  parameter 
as  the  adjustable  parameter  which  is  used  to  force  MINIMUF  to  give  the 
measured  MOF  value.  This  effective  10.7  cm  flux,  wnich  we  shall  call  simply 
the  "effective  flux",  is  usually  different  from  direct  measurements  of  the 
10.7  cm  flux,  ana  therefore  loses  its  literal  meaning.  It  can  nevertheless  be 
considered  as  an  "adjusted  10.7  cm  flux"  which  is  (in  principle)  independent 
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or  geographical  position  and  therefore  suitable  as  a  parameter  with  which  the 
correlation  of  ionospheric  properties  may  be  tested. 

This  technique  ot  correlating  an  adjustable  moael  parameter  is  somewhat 
diflerent  trom  previous  methods  in  wnich  differences  between  model  values'and 
measured  values  ot  MuF  were  correlated  Ouisn  J.97b;  Gautier  ana  Lacharisen 
19b5;  beckwitn  ana  Kao  1975). 

Correlation  of  effective  Fluxes  on  a  rair  of  Faths 

To  illustrate  the  correlation  ot  effective  fluxes  for  a  pair  ot  paths  we 
use  data  from  the  Solid  Snield  experiment  Ulffelman  and  Harnish  1982),  a  test 
of  HF  communications  procedures  which  occurred  on  3-20  May  1981.  These  data 
were  chosen  because  ot  the  relatively  large  quantity  of  data  availaole:  there 
are  lb04  MOF  measurements,  covering  401  hours  of  sounder  operation  spread  over 
17  days.  MOFs  were  measured  on  the  six  paths  listed  in  Table  1  and 
illustrated  in  Figure  1^ 

Figure  2  shows  a  comparison  of  tne  eftective  fluxes  measured  for  the 
MacDill-Norfoln  and  Le jeune-Norfolk  paths.  Each  of  the  three  plots  includes 
data  acquired  during  a  specific  nour  of  each  of  a  number  of  days.  in  Figure 
2a,  for  example,  the  nine  data  points  refer  to  data  acquired  between  0700  and 
ObOOZ  on  each  of  nine  days.  Since  ionograms  were  made  every  fifteen  minutes, 
each  effective  flux  value  plotted  is  tne  flux  for  which  MINIMUF  produces  MUF 
values  closest  to  the  measured  values  for  the  tour  measurements  of  the  nour, 
in  the  sense  of  smallest  rms  difference.  Effective  fluxes  were  calculated 
with  the  NKL  OISI  computer  program.  Figures  2a,  2b,  and  2c  were  chosen  to 
illustrate  situations  in  which  the  correlation  is,  respectively,  about 
average,  a  little  better  tnan  in  most,  and  worse  than  in  most  situations. 

For  each  of  the  plots  in  Figure  2  a  linear  regression  line  in  the  form 

Y  -  a  ♦  b  X.  (1) 

was  calculated,  where  X  and  Y  measure  the  effective  fluxes  on  Che 
Le jeune-Norfolk  and  MacDill-horfolk  patns,  respectively.  The  fitted  line  may 
be  used  with  M1N1MUF  for  forecasting,  since  a  measured  MOF  and  the  derived 
effective  flux  on  the  Le jeune-Norfolk  path  lead,  through  equation  1,  to  an 
effective  flux  and  then  to  a  MINIMUF  MOF  prediction  on  the  MacDill-Norfolk 
patn.  This  technique  of  using  a  model  and  a  nut'  measurement  on  one  path  to 
forecast  a  MOF  on  another  path  has  been  termed  "Updating"  and  has  been  the 
subject  of  a  number  of  NkL  studies  lUffelman  1981;  Utfelman  et  al  1981,  1982, 
1984a,  1984b).  Me  note  that  the  use  of  equation  l  to  forecast  the  effective 
flux  on  an  unmeasured  path  represents  a  departure  from  previous  applications 
of  Updating,  where  the  simpler  relation  Y  ■  X  was  used.  In  principle,  Y  could 
be  calculated  as  a  polynomial  in  X  of  arbitrarily  high  degree,  but  the 
presently-available  accuracy  of  MOF  and  effective-flux  measurements  doesn't 
justify  the  nigher  precision. 

In  some  correlation  plots  it  was  noticed  that  one  or  two  "wild"  points 
resulted  in  a  linear  regression  line  with  a  slope  differing  greatly  from  the 
expected  slope  of —45°,  as  illustrated  in  Figure  3.  While  one  must  be 
cautious  about  throwing  out  data,  occasional  events  such  as  violent 
ionospneric  disturbances,  temporary  electrical  disturbances,  or  blunders  in 
measurement  can  sometimes  obscure  a  real  correlation.  Therefore,  we  have  in 
some  cases  processed  the  data  twice,  lirst  usin6  all  the  data,  and  then 
eliminating  by  subjective  exclusion  up  to  three  points  per  plot.  Not 
surprisingly,  the  exclusion  process  can  greatly  improve  tne  correlation  plot, 
as  illustrated  in  Figure  3.  With  a  larger  data  base  this  procedure  would 
presumably  become  unnecessary. 
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The  correlation  of  erfective  fluxes  on  a  pair  or  patns  can  be  measured 
quantitatively  Dy  tne  correlation  coefficient  r.  Let  ,  X2,  ...  Xr  be 
n  values  or  effective  llux  over  path  1  at  tne  times  t^,  t2,  ...  tn.  Let 

^li  ^2*  •♦•^n  cne  values  of  effective  flux  over  path  Z  at  the  same 

times.  The  correlation  of  the  effective  fluxes  on  the  two  paths  is  defined  by 

the  correlation  coefficient  r: 


n  n  n 


(.See,  for  example,  Kef.  8,  p.  523). 

Thus  defined,  the  correlation  coefficient  measures  the  degree  of  linear 
relationship  between  the  variables  X  and  Y  and  is  independent  of  the  units  in 
wnich  the  variables  are  measured.  It  can  vary  between  ral  (.perfect 
correlation)  to  -1  (perfect  anti-correlation}.  When  applied  to  forecasting, 
ral  implies  forecasting  with  near-certainty,  ana  r"o  implies  uselessness.  A 
negative  correlation  coefficient  is  inconsistent  witn  our  premises;  it  is 
considered  a  non-physical  result  anu  therefore  in  error. 

The  significance  of  the  correlation  coefficient  must  be  measured  in  terms 
of  its  ability  to  forecast  the  value  of  effective  flux  on  one  pach  when  the 
effective  flux  on  the  other  path  has  been  measured.  To  make  such  inferences, 
one  must  assume  tnat  vanaoles  X  and  Y  are  both  normally  distributea.  buch  an 
assumption  seems  reasonable,  though  we  don't  verify  that  here.  A  useful 
interpretation  is  that 

r^  •  Variation  in  Y  explained  by  regression  line. 

Total  variation  observed  in  Y. 

The  quantity  r^  is  xnown  as  the  coefficient  of  determination.  Thus 
correlation  coefficients  on  the  oraer  of  U.9  or  greater  imply  r2>0.8,  which 
indicates  that  the  assumption  of  linear  relationship  between  the  effective 
fluxes  is  reasonably  valid. 

Time  Variation  of  Correlation  Coefficient 

We  have  computed  the  correlation  coefficient  as  a  tunction  of  time  tor 
eighteen  combinations  of  the  six  paths  involved  in  the  Solid  Shield  project. 
The  results  are  shown  in  Figure  4.  For  most  pairs  of  paths  tnere  are  two 
plots:  one  uses  all  the  data,  and  the  other  uses  data  from  which  "wild 
points"  have  been  discarded.  We  make  tne  following  observations: 

1.  There  are  many  measurements  of  correlation  coefficients  between  0.9 
and  1,  corresponding  to  coerficients  of  determination  between  0.8  ana  1,  which 
implies  a  reasonable  likelihood  of  using  the  correlation  of  effective  fluxes 
as  a  predictive  tool; 

2.  In  some  cases,  however  (such  as  in  the  Le jeune-uriver ,  uriver-Norf oik 
correlation,  plot  2),  the  improvement  in  correlation  caused  by  subjectively 
eliminating  "wild  points"  is  striking.  Tne  eliminated  data  in  this  situation 
were  a  string  of  consecutive  data  points  on  a  single  day,  probably  indicating 
a  systematic  error  in  measuring  tne  M OF  on  one  of  the  paths  during  that  day. 
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3.  Seven  of  the  diurnal  correlation  plots  (.Figure  4,  plots  1,  8-13J 
indicate  a  large  decrease  in  tne  correlation  coefficient  curing  a  period  of 
four  nours  centered  on  approximately  lbOO  liT.  Since  all  involve  the 
Macui 1 i-Nor folk  patn,  tne  cause  is  probably  some  systematic  error  in  tnose 
data . 

4.  t'siae  from  the  gross  features  mentioned  above,  the  correlation 
coefficient  tends  to  be  at  its  nighest  around  midnight  UT  (1900  LT)  and  tends 
to  be  jomewhat  oecreaseu  around  1400  OT  and  possibly  around  0300  01  (0900  ana 
1200  L.n,  respectively).  These  trends  are  not  well  defined,  however. 

Diurnal  Variation  of  Regression  Coefficients 

figure  3  shows  the  diurnal  variation  of  the  constants  A  and  B  determined 
for  the  linear  regression  fit  to  the  effective  fluxes  calculated  for  the 
Le jeune-Norfolk,  Driver-Bragg  pair  of  paths.  The  hand-drawn  lines  show  that 
the  coefficients  may  possibly  be  oescribed  as  varying  roughly  sinusoidally, 
with  a  twelve-hour  period.  Sucn  a  regularity  would  tend  to  indicate  a 
deficiency  in  MINIMUF's  ability  to  account  ror  geographical  MuF  dependence,  at 
least  as  far  as  this  particular  example  is  concerned. 

Similar  plots  for  several  otner  pairs  of  patns  were  maae  to  see  if  tne 
regularity  observed  for  the  Le jeune-Norfolk,  Driver-Bragg  regression 
coefficients  is  typical.  Tne  resulcs,  shown  in  figure  b,  are  ambiguous.  From 
tnese  data,  one  can  only  conclude  that  in  some  instances  a  periodic  variation 
in  regression  coetficients  occurs,  and  that  a  model  improvement  should  be 
possible. 

Spatial  De-correlation 

It  is  reasonable  to  expect  the  correlation  of  effective  fluxes  determined 
for  a  pair  of  paths  to  decrease  as  the  distance  between  the  path  centers 
(called  control  points)  increases.  To  test  tnia  effect,  tne  average 
correlation  coefficient  for  each  of  the  eighteen  pairs  of  paths  used  in 
creating  Figure  <♦  is  plotted  as  a  runction  ot  separation  ot  tne  control  points 
(Figure  7).  as  expected,  the  correlation  decreases  with  increasing 
control-point  separation.  At  longer  path  separations  the  scatter  in  the 
correlation  coefficient  becomes  large.  Because  of  the  limited  scope  of  the 
Solid  bhiela  test  (there  are  no  control-point  separations  greater  than  700 
km),  it  is  not  possiole  to  confirm  tne  decrease  in  f0F2  correlation 
measured  by  Kush  (197b)  for  ionospheric  position  separations  up  to  5000  km. 

ill.  mof  Furl casts 

We  consider  three  ways  to  forecast  the  MUF,  all  of  them  using  MINIMUF: 

1.  MINIMUF  alone 

2.  HINIMUF  updated  by  sounder  measurements  on  the  same  path 

3.  MIN LMliF  updated  by  sounder  measurements  on  a  different  path. 

forecasting  by  MInlMUF  Alone 

Given  a  specified  communications  path,  day  of  year,  and  solar  lu.7  cm  flux 
value,  the  MINIMUF  model  predicts  the  MOF  as  a  funcion  of  time  of  day.  There 
is  some  oeoate  about  the  most  appropriate  value  of  solar  flux  to  use;  in  some 
cases  tne  daily  value  (as  broadcast  hourly  by  WWV)  seems  to  work  well,  while 
in  other  cases  a  value  equal  to  the  flux  averaged  over  the  previous  five  days 
is  preferred.  To  continue  with  the  previously-discussed  example  from  the 
bolio  bnieid  test  (Le jeune-Nor folx  path,  3-19  May  1981),  we  have  calculated 
the  MINIMUF  MOF  predictions  and  compared  them  with  the  observed  MUFs.  These 
data  include  approximately  350  nours  of  measurements  spanning  a  period  or 
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seventeen  days.  The  rms  error  of  tnese  forecasts  turned  out  to  be  2.49  MHz 
when  tne  daily  soiar  ilux  value  wa9  useu,  ana  2. bn  Mnz  wnen  the  rive-aay 
average  was  used. 

These  values  are  inciuaed  in  figure  8,  wnicn  shows  the  rms  difierence 
between  the  forecast  and  measured  MOF  values  as  a  function  of  the  time  delay 
between  the  time  of  forecast  ana  the  time  of  measurement,  bince  the  MINIMUF 
forecast  changes  only  once  per  day,  when  a  new  solar  flux  value  becomes 
available,  the  rms  error  is  independent  of  tne  delay  between  tne  times  of 
forecasting  and  measurement,  for  times  up  to  twenty-four  hours.  Thus,  the 
MlhlMbr -alone  rms  error  appears  as  a  norizontal  line  in  figure  8. 

Forecasting  by  the  Updated  MINIMUF  Model 

The  prospect  of  improving  torecasting  accuracy  by  "Updating"  is  basea  on 
tne  presumption  of  low-trequency  (f£l  day  variations  in  tne  MoF  which 
are  not  accounted  for  in  the  MINIMUF  model.  If  Che  model  MOf  forecast  is 
multiplied  by  a  constanc  wnicn  maxes  the  result  agree  perfectly  with  a  current 
measured  moF,  it  is  reasonable  to  expect  that  multiplying  Che  model  forecast 
for  a  later  time  by  the  same  constant  will  improve  the  forecast,  bince  the 
solar  10.7  cm  flux  enters  into  the  MINIMUF  calculations  essentially  as  a 
multiplicative  factor,  it  has  been  convenient  to  modify  the  model's 
predictions  oy  using  an  "effective  flux"  instead  of  the  actual  one.  The 
etfeccive  flux  is  determined  at  the  time  of  Updating  and  is  used  to  obtain 
MliilMUF  forecasts  until  a  new  effective  flux  is  determined.  The  purpose  of 
this  section  is  to  snow  an  example  of  Updating  methods  and  to  show  how  its 
efficacy  in  improving  MOF  forecasting  can  be  measured. 

The  eftective  flux  may  be  determined  from  sounder  measurements,  such  as 
the  ionograms  from  the  oblique-incidence  souuders  used  in  the  bolid  Shield 
test,  nfrective  flux  measurements  may  De  made  either  on  the  path  for  which 
forecasts  are  desired  or  on  some  other  path.  When  using  an  effective  flux 
measured  on  a  path  ditferent  from  tne  patn  on  which  tne  forecast  is  desired, 
it  is  necessary  to  use  a  correlation  function  similar  to  those  illustrated  in 
Figure  2  to  determine  the  effective  flux  on  the  path  tor  wnich  predictions  are 
being  made. 

Figure  8  shows  the  rms  error  between  forecast  and  measured  rtuF.-i  on  the 
Lcjeune-Norfolk  path  as  a  function  of  the  time  delay  between  forecast  and 
measurement.  One  curve  results  from  predictions  using  effective  flux 
measurements  from  the  same  path.  The  calculations  for  this  curve  include  all 
data  for  that  path  for  which  botn  torecast  and  measurement  (.within  24  hours ) 
were  possible.  In  this  case,  the  correlation  measurements  are  not  involved. 
The  otner  two  curves  are  derived  from  forecasts  based  on  effective  fluxes 
measured  on  two  other  paths:  Driver-Ft.  dragg,  and  MacDill-Ft.  Bragg.  In 
these  cases,  the  corresponding  eftective  fluxes  for  the  Le jeune-Nortolk  patn 
were  determined  from  correlation  plots  of  the  type  illustrated  in  Figure  2. 

Tne  presentation  of  data  in  Figure  6  anu  conclusion  cased  on  the  figure  are 
flawed  by  the  fact  that  the  MuF  measurements  used  to  check  the  forecasts  were 
also  inciuaed  in  the  aaca  useu  to  determine  the  correlation  coet f icients . 

With  a  larger  data  base,  it  would  be  possible  to  avoid  this  difficulty. 

Figure  8  is  included,  however,  because  it  typifies  the  data  required  to 
determine  the  efficacy  of  Updating  techniques. 

Comments  on  Updating  Techniques 

The  following  observations  can  oe  made  on  the  previous  example  ot  Updating 
procedures : 


6 


a 


1.  For  che  Lejeune-Nortolk  aata  presenced  in  Figure  d,  the 
Updated-MINIMUF  forecasts  for  delay  times  up  to  about  six  hours  result  in  a 
smaller  rms  error  chan  forecasts  using  MIaL*lUF  alone. 

2.  For  prediction  times  between  twelve  and  eighteen  hours,  Figure  8  shows 
tne  somewnat  surprising  result  chac  tne  Updated-MINlinUF  forecasts  are  interior 
co  chose  using  MINLMUF  alone. 

3.  There  is  a  martced  similarity  in  the  Updatec-MlM LFiUF  results  tor  tne 
chree  ways  of  determining  Che  effective  flux  on  the  Le jeune-Norfoik  path.  it 
seems  surprising  tnat,  except  for  snore  delay  times  (,1-b  nours)  between 
forecast  ana  observation,  Che  forecasts  based  on  Le jeune~Norf oik  measurements 
are  interior  to  chose  basea  on  measurements  trom  che  ocher  paths. 

k.  The  MINIMUF-alone  rms  error  of~2.6  MHz  are  substantially  lower  than 
Che  generaily-applicable  3.3  MHz  errors  claimed  by  rtlNIriUF.  It  is  possible 
chat  che  Lejeune -Nor folk  path  is  a  relatively  favorable  case  for  the 
application  of  MINIMUF  prediction. 

3.  There  is  a  trend  for  forecasts  tor  '■>•12  hours  ahead  to  be  worse  than 
those  for  ~  2k  hours  ahead,  regardless  of  the  type  of  Updating.  This  is 
prooably  inherent  in  che  moael;  use  of  a  more  sophisticated  model  might 
eliminate  this  effect. 

6.  It  should  oe  notea  that  the  results  presented  in  this  report  are  cased 
on  a  very  small  data  sample  and  may  not  represent  general  ionospheric 
conditions.  However,  the  procedures  used  are  those  which  may  be  used  to 
assess  the  validity  of  Updating  Cor  other)  techniques  of  MOF  forecasting. 

7.  The  dependence  of  forecast  accuracy  as  a  function  of  time  of  torecast 
has  not  been  explored.  For  example,  it  may  be  that  forecasts  made  at  sunrise 
for  i2  ahead  are  oetcer  chan  these  made  at  noon  for  12  hours  aheaa.  The  boiid 
Shield  data  seemed  insufficent  in  quantity  for  an  adequate  discussion  of  this 
dependence 

Although  tne  present  work  has  used  a  single  parameter  Uhe  effective  flux) 
measured  on  a  sounded  path  to  update  the  model  forecast  of  the  same  parameter 
on  tne  same  patn  or  on  an  unsounded  path,  more  elaborate  schemes  are 
possible.  For  example,  the  model  forecast  of  a  parameter  may  be  updated  using 
measurements  of  che  same  parameter  at  a  number  ot  positions  and  times,  or  it 
may  be  updated  with  recent  measurements  of  a  numoer  of  different  parameters 
lUautier  and  Zacnansen  1963,  Zacharisen  i9o3,  Zacharisen  ana  Gautier  l9t>4). 
Updating  with  a  large  quantity  of  sounder  data  would  be  appropriate  for  a 
network  of  sounaers  established  for  forecasting  propagation  conditions  over  a 
large  geographical  region. 

IV.  FRfcCjlitNCY  MANAGER'S  e OX-FORECASTING  PKOCEDUKt 

Management  of  an  HF  communications  network  demands  the  assignment  of 
frequencies  in  such  a  way  tnat  the  requirements  ot  connectivity,  data  races, 
covertness,  etc.,  are  achieved.  Here  we  are  concerned  witn  only  one  of  the 
constraints  to  che  assignment  process:  the  requirement  that  assigned 
frequencies  be  at  or  below  the  MOF  for  the  assigned  communications  path. 

Given  a  measurement  of  MoF  on  a  sounded  patn,  we  attempt  to  produce  a 
procedure  for  determining,  with  a  predetermined  probability  of  success,  the 
highest  frequency  which  does  not  exceed  che  MOF  on  the  communications  path. 
Such  a  frequency  may  be  called  the  Frequency  of  Optimum  Transmission,  or  FOX. 
It  is  not  a  single  quantity,  since  it  depends  on  the  specified  probability. 

We  thus  speak,  for  example,  of  the  "902  FOT"  for  a  given  path  and  time, 
meaning  tne  highest  rrequency  which,  witn  902  probability,  will  be  below  the 
actual  MOt .  This  definition  differs  from  that  used  in  some  frequency 
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management  procedures,  wnere  Che  MUF  is  defined  as  a  percentage  (typically 
oo*,)  of  the  MOF  (uavies  1965).  Utners  simply  define  tne  FOT  as  the  90X  FOT 
aetinea  aoove  (uautier  and  Zacnansen  i9ho>. 

Probability  Interval  of  an  Effective  Flux  Forecast 

In  a  previous  seccion  it  was  snown  now  tne  correlation  of  effective  tluxes 
on  two  patns  could  De  approximated  by  a  linear  regression  line  and  measured  by 
tne  coerficient  of  correlation.  lhen,  in  discussing  Upuating,  tne  measured 
effective  flux  on  a  reference  path  and  the  known  correlation  information  were 
used  co  rorecasc  tne  effective  tiux  on  a  communication  patn.  Here  we  show  the 
degree  of  confidence  with  which  this  forecast  may  be  used,  given  a  specified 
amount  of  correlation  data. 

We  assume  that  a  series  of  n  effective  fluxes  X^,  X^»  •••  X^  on  one 
patn  (.called  the  control  patn,)  ana  tne  corresponding  efrective  fluxes  on  F^, 
F2,  ...  rn  on  a  second  patn  (called  tne  communications  path)  nave  been 
measured  ana  that  a  linear  regression  line  to  approximate  the  correlation 
lines  uas  been  calculated.  Thus,  given  a  subsequent  ettective  flux  X  measured 
on  the  control  path,  the  best  estimate  of  effective  flux  on  the  communications 
path  is 

A 

r  ■  a  0  X 

Inis  is  not  an  appropriate  efrective  flux  to  use  for  frequency  management 
purposes  because  there  is  a  S0/«  probability  that  the  FOT  calculated  with  this 
erfective  flux  will  exceed  tne  actual  M<JF.  satisfactory  communications 
require  a  higher  probability  of  success  -  perhaps  80  to  95X. 

The  prooable  error  in  forecasting  an  effective  flux  rrom  the  measured 
correlation  relationship  may  be  calculated  by  appropriate  statistical 
cecnniques  (see  Harnett,  chapter  10  for  a  goou  discussion).  For  values  of 
control  patn  effective  flux  near  tne  mean  value 

n 

*  ■  l 

n  1.1 

the  standard  error  of  forecast  in  the  corresponding  values  or  F  is 


A 

wnere  rj_  a  a  +  b  X^  is  the  linear  regression  estimate  of  F  corresponding 
to  tne  control  path  enective  ilux  XL.  ror  values  of  control  path  erfective 
flux  X  not  near  the  mean  the  standard  error  of  forecast  is  modified  as  follows 


The  quantity  under  cne  radical  sign  quantifies  tne  greater  inaccuracy 
associated  witn  predictions  of  F  for  control  path  effective  flux  values  X  far 
from  the  mean. 
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A 

"hue  cne  Linear-regression  estimate  F  tor  cne  effective  flux  gives  the 
best  prediccion  possible  from  cne  n  available  measurements,  we  can  also 
calcuiace  cne  likeiinood  cnat  cne  actual  value  will  depart  trom  Che  estimated 
value  by  a  specified  amount.  Estimates  in  the  present  case  (i.e.,  a  case  in 
whicn  Cne  standard  deviation  for  all  possible  values  of  r  is  unknown!  are 
expected  to  follow  a  t-distribution  with  v  a  n-2  degrees  of  treedom.  That  is, 
Che  prooabiiity  tnat  the  actual  effective  flux  will  depart  from  the  estimated 
value  F  by  the  amount  t  is  given  by  P  It,  v),  a  function  which  can  be  eitner 
computed  or  reau  from  a  table.  (.Most  often,  the  integrals  of  P  over  specific 
ranges  of  t  are  tabulated.)  The  normalized  variable  t  is  the  effective  flux 
measured  in  units  of  tne  standard  error  of  forecast  Sf.  Thus  a  04 
confidence  interval  for  Che  actual  flux  may  be  defined  by  the  extremes 

F  1  CC,VSf, 

where  tcv  is  cne  numoer  tor  whicn 


+tc ,  v 

P  ( t , v)  dt  =  C/100 

,v 

Figure  9  illustrates  cne  application  of  this  procedure  to  the 
correlation  data  collected  on  eleven  days  in  the  time  interval  0300  -0400  2  on 
tne  two  Solid  bhield  paths  MacDill-Norfolk  and  Le jeune-Nor tolk.  The 
correlation  between  the  effective  flux  F  on  the  MacDill-Norfolk  path  to  the 
effective  flux  X  on  the  Le jeune-Norfolk  patn  is  expressed  by  the  linear 
regression  line 
A 

F  [Mhz)  =  2d. 8b  *  U.83L  X  iMKzj. 


The  80,  90  and  95  percent  confidence  intervals  calculated  by  the  above 
equations  are  also  plotted.  Inus,  tor  example,  a  measured  eftective  flux  on 
tne  control  patn  (Lejeune-Norfolk)  results  in  a  best  estimate  for  che 
effective  tlux  on  the  MacOil 1-norf oik  patn  of  137.0,  with  an  804  probability 
that  the  actual  effective  flux  will  be  between  118  and  57  and  a  954 
probaoilicy  that  it  will  fall  oecween  10b  and  lo9. 

FOT  Forecasting  Procedure 

The  MIuIMOF  program  may  be  used  to  transform  tne  forecasts  of  effective 
flux  and  confidence  intervals  on  the  MacDill-Norfolk  patn  (.Figure  9)  into  FOT 
forecasts,  thereby  maxing  the  data  more  readily  usable  tor 

frequency-management  purposes  (see  Figure  10).  In  transforming  the  effective 
flux  forecasts,  only  tne  linear  regression  line  and  the  Lower  limits  of  the 
probability  intervals  are  retained  and  their  labels  are  changed  to  correspond 
to  a  new  interpretation.  The  new  curves  represent  the  most  probable  value  for 
the  MOF ,  and  the  frequencies  which  are  expected  to  be  below  the  actual  MOF 
with  probabilities  of  904  954,  ana  9  7 . 54. .  The  latter  curves  are  thus  labeled 
as  "FOTs"  (Frequencies  of  Optimum  Transmission),  and  the  most  probable  MOF 
curve  is  now  che  "504  FuT".  This  set  ot  information  snouid  be  of  more 
assistance  in  selecting  a  propagating  frequency  than  a  single  value  of  MUF 
whose  relation  to  the  variance  in  the  expected  MOF  is  unknown. 
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Table  1  —  Solid  Shield  Propagation  Paths 


it 

1 

TO  FT . 3RAGG , 

NC 

LATITUDE 

+35.15 

LONGITUDE 

-78.98 

FROM 

LATITUDE 

LONGITUDE 

DISTANCE (km) 

MIDPOINT 

LAT .  LONG . 

* 

DRIVER,  VA 

+36.82 

-76.50 

290 

+35.99 

-77.75 

HURLRERT ,  FL 

+30,3 

-86.4 

879 

+32.78 

-82.79 

SHAW  AFB ,  SC 

+34.97 

-80.48 

138 

+35.06 

-79.73 

1 

MACDILL  AFB, 

FL 

+27.85 

-82.48 

877 

+31.51 

-80.80 

TO  NORFOLK,  VA 

LATITUDE 

+36.67 

LONGITUDE 

-76.23 

1 

FROM 

LATITUDE 

LONGITUDE 

DISTANCE (km) 

MIDPOINT 

LAT .  LONG . 

* 

DRIVER,  VA 

+36.82 

-76.50 

29 

+36.75 

-76.36 

HURLBERT,  FL 

+30.3 

-86.4 

3.178 

+33.59 

-81.50  . 

\ 

SHAW  AFB,  SC 

+34.97 

-80.48 

427 

+35.84 

-78.38 

MACDILL  AFB,  FL 

+27.85 

-82.48 

1143 

+32.30 

-79.51 

LEJEUNE ,  NC 

+34.67 

-77.35 

244 

+35.67 

-76.80 

EFFECTIVE  10.7  CM  FLUX  (MscDM  to  Norfofc) 
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Fig.  2  -  Three  examples  of  correlation 
of  effective  fluxes  on  the  MacDill- 
Norfolk  and  LeJeune-Norfolk  paths.  An 
’effective  flux”  is  the  10.7  cm  solar  flux 
for  which  the  MINIMUF  MUF  predic¬ 
tion  equals  an  observed  MUF.  Examples 
a,  b,  and  c  contain  data  for  the  1-hour 
intervals  starting  0800Z,  0300Z,  and 
1600Z,  respectively,  within  the  Solid 
Shield  test  period  8  •  18  May  1981. 
Straight  lines  are  least-squares  regression 
fits  to  the  data. 


EFFECTIVE  FLUX  (LoJeun#  to  Norfolk) 


Fig.  3  —  Effect  of  removing  "wild  points"  from  the  correla¬ 
tion  data.  Both  plots  show  the  correlation  of  effective 
fluxes  on  the  MacDill-Norfolk,  LeJeune-Norfoik  paths  for 
the  hour  1500-1600Z.  In  the  lower  plot,  two  "wild  points" 
have  been  subjectively  removed.  Straight  lines  are  least- 
squares  regression  fits  to  the  data. 


Fig.  4  —  Correlation  coefficient,  as  a  function  of  diurnal  time, 
measuring  the  correlation  of  effective  fluxes  for  each  of  eigh¬ 
teen  pairs  of  Solid  Shield  propagation  paths.  For  each  pair,  the 
plot  on  the  left  was  produced  using  all  the  available  data;  in  the 
plot  on  the  right,  some  data  points  were  calculated  after  subjec¬ 
tively  deleting  'wild  points”  from  the  data  (see  text). 


Fig.  4  (Continued)  —  Correlation  coefficient,  as  a  function  of 
diurnal  time,  measuring  the  correlation  of  effective  fluxes  for 
each  of  eighteen  pairs  of  Solid  Shield  propagation  paths.  For 
each  pair,  the  plot  on  the  left  was  produced  using  all  the  avail' 
able  data;  in  the  plot  on  the  right,  some  data  points  were  calcu¬ 
lated  after  subjectively  deleting  "wild  points”  from  the  data  (see 
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Fig.  4  (Continued)  —  Correlation  coefficient,  as  a  function  of 
diurnal  time,  measuring  the  correlation  of  effective  fluxes  for 
each  of  eighteen  pairs  of  Solid  Shield  propagation  paths.  For 
each  pair,  the  plot  on  the  left  was  produced  using  ail  the  avail* 
able  data;  in  the  plot  on  the  right,  some  data  points  were  calcu¬ 
lated  after  subjectively  deleting  "wild  points”  from  the  data  (see 
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Fig.  5  -  Diurnal  variation  of  the  linear  regression  coefficients  describing  the  correlation  of 
effective  flux  F  on  the  driver  to  ft.  Bragg  path  with  the  effective  flux  X  on  the  Lejeune  to  Nor¬ 
folk  path.  The  linear  regression  lines  is  F  -  A  +  BX. 
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Fig.  7  —  Decrease  in  correlation  coefficient  with  distance  between  the  control 
points  for  18  pairs  of  Solid  shield  paths.  The  control  point  is  the  geographical  mid¬ 
point  of  a  communication  path.  (It  is  assumed  that  all  propagation  is  by  single-hop 


i — i — i — i — i — i — i — i  i  i  i  i  i  i  i  i  i  i  ;  i  i  i  r 


PREDICTION-MEASUREMENT  DELAY  (HOURS) 

Fig.  8  —  RMS  error  between  actual  MOFs  on  the  Lejeune  to  Norfolk  path  and  MOFs  predicted 
by  a  variety  of  techniques,  as  a  function  of  time  delay  between  prediction  and  measurement. 
Data  include  approximately  335  hourly  predictions  for  each  value  of  time  delay. 
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EFFECTIVE  FLUX  ON  LEJEUNE-NORFOLK  PATH 

Fig.  9  —  Confidence  intervals  in  the  prediction  of  effective  flux  on  the  MacDill  to  Norfolk  path 
in  the  time  interval  0300Z-0400Z,  using  MIN1MUF  updated  with  Lejeune  to  Norfolk  data. 
Given  a  measurement  of  effective  flux  on  the  Lejeune  to  Norfolk  path,  the  best  estimate  of 
effective  flux  on  the  MacDill  to  Norfolk  path  is  given  by  the  linear  regression  line.  Three  pairs 
of  curves  represent  the  intervals  within  which  the  MacDill  to  Norfolk  effective  flux  may  be 
expected  to  fall  with  probabilities  of  30,  90,  and  95  percent. 
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EFFECTIVE  FLUX  ON  LEJEUNE-NORFOLK  PATH 

Fig.  10  —  Curves  for  forecasting  the  Maximum  Usable  Frequency  (MUF)  on  the  MacOill  to 
Norfolk  path  for  the  hour  0300Z-0400Z.  Given  a  measurement  of  effective  flux  on  the 
Lejeune  to  Norfolk  path,  the  upper  curve  gives  the  best  estimate  of  the  MOF  on  the  MacDill  to 
Norfolk  path.  The  other  three  curves  are  guides  for  determining  a  MUF  with  a  specified 
probability  that  this  frequency  will  be  at  or  below  the  actual  MOF.  For  example,  the  "80% 
MUF  is  the  highest  frequency  which,  with  80%  probability,  will  be  at  or  below  the  actual  MOF. 
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